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Abstract 

First evidence of the B*^ — > J/4) cu decay is found and the — ;> 3/^\)r\ and B^ — > 3/^\)r\' 
decays are studied using a dataset corresponding to an integrated luminosity of 
1.0 fb~^ coUected by the LHCb experiment in proton-proton cohisions at a centre-of- 
mass energy of y/s = 7 TeV. The branching fractions of these decays are measured 
relative to that of the B*^ — )• J/i}; p*^ decay: 

|g^= 0.89 ±0.19 (Stat) iO:?^(syst), 

|g^= 14.0 ± 1.2 (Stat) t\i (syst) 

where the last uncertainty is due to the knowledge of /d//s, the ratio of b-quark 
hadronization factors that accounts for the different production rate of B'^ and Bg 
mesons. The ratio of the branching fractions of Bg — )• J/4'Tl' E^iid Bg — )• J/il^rj decays 
is measured to be 

S(B0 ^ J/t|;ti') _ , +0.06 



S(BO ^ J/^\>ri) 



0.90 ± 0.09 (stat) (syst). 
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1 Introduction 



Decays of B mesons into a J/i|) and a light meson are dominated by color-suppressed tree 
diagrams involving b — )■ ccs and b — ccd transitions (see Fig. Contributions from 
other diagrams are expected to be small p]. Measurements of the branching fractions of 
these decays can help to shed light on hadronic interactions. The decay B° — )■ J/ip tu has 
not been observed previously. The CLEO collaboration has set the most restrictive upper 
limit to date of i3(B° J/^\)^v) < 2.7 x 10"^ at 90% confidence level g. 

The B° ^ J/i|;riW decays were observed by the Belle collaboration [3] with branching 
fractions i3(B° ^ J/ipii) = (5.10 ± 0.50 ± 0.25 IJ;^^) x 10"^ and i3(B° ^ J/\ljri') = 
(3.71 ±0.61 ±0.18 

-0^7) ^ 10 ^) where the first uncertainty is statistical, the second is 
systematic and the third one is due to an uncertainty of the number of produced B^B^ 
pairs. Since both final states are CP eigenstates, time-dependent CP violation studies and 
access to the Bg — B^ mixing phase cjjs will be possible in the future |4j. The theoretical 
prediction for these branching fractions and their ratio relies on knowledge of the t) — rj' 
mixing phase cf)?. Taking = (41.4 ±0.5)° [5j and ignoring a possible gluonic component 
and corrections due to form factors, the ratio becomes 

— 1 Ofi +0.10 

S(BO ^ J/^\>^) " tan2(t)p " ^'^^-o-os- 

Here J^' ' is the phase space factor of the B^ J/i|) r\^'^ decay and the uncertainty is due 
to the inaccuracy in the knowledge of the mixing phase. As discussed in Ref. [l], a precise 
measurement of this ratio tests SU{3) flavour symmetry. In addition, in combination with 
other measurements, the fraction of the gluonic component in the r\' meson can eventually 
be estimated [6|. 



b c b c 




d d s s 



Figure 1: Examples of the dominant diagrams for the B^^^ — )■ J/i|)X° decays (where 
X° = ri, ri', tu or p°). 

The analysis presented here is based on a data sample corresponding to an integrated 
luminosity of 1.0 fb~^ collected by the LHCb detector in 2011 in pp collisions at a centre- 
of-mass energy of ^/s = 7 TeV. The branching fractions of these decays are measured 

relative to B{B^ — )• J/i|) p°) and the ratio jg^j^o _^ determined. 
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2 LHCb detector 



The LHCb detector [7| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r/ < 5, designed for the study of b- and c-hadrons. The detector includes a high 
precision tracking system consisting of a silicon-strip vertex detector surrounding the pp 
interaction region, a large-area silicon-strip detector located upstream of a dipole magnet 
with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw 
drift tubes placed downstream. The combined tracking system has a momentum resolution 
Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, and an impact parameter 
resolution of 20 |am for tracks with high transverse momentum (pt)- Charged hadrons 
are identified using two ring-imaging Cherenkov (RICH) detectors. Photon, electron and 
hadron candidates are identified by a calorimeter system consisting of scintillating-pad and 
pre-shower detectors, and electromagnetic and hadron calorimeters. Muons are identified 
by a system composed of alternating layers of iron and multiwire proportional chambers. 

The trigger consists of a hardware stage, based on information from the calorimeter and 
muon systems, followed by a software stage which applies a full event reconstruction. This 
analysis uses events triggered by one or two muon candidates. In the case of one muon, 
the hardware level requirement was for its to be larger than 1.5 GeV/c; in case of two 
muons the restriction ^/Pti ■ Pt2 > 1-3 GeV/c was applied. At the software level, the two 
muons were required to have an invariant mass in the interval 2.97 < m^+^- < 3.21 GeV/c^ 
and to be consistent with originating from the same vertex. To avoid the possibility that 
a few events with high occupancy dominate the trigger processing time, a set of global 
event selection requirements based on hit multiplicities was applied. 

For the simulation, pp collisions are generated using Pythia 6.4 [s] with a specific 
LHCb configuration \9\. Decays of hadronic particles are described by EvtGen jlO] in 



which final state radiation is generated using Photos 11 . The interaction of the generated 



particles with the detector and its response are implemented using the Geant4 toolkit [12J 



as described in Ref. 13 . The digitized output is passed through a full simulation of both 



the hardware and software trigger and then reconstructed in the same way as the data. 



3 Data sample and common selection requirements 

The decays B^^^ J/v|jX° (where X° = r|, r[', w and tt+tt^) are reconstructed using 
the J/4) — )■ \Ji^\Ji~ decay mode. The X° candidates are reconstructed in the r\ — >■ yy, 
r\ —7- 7t+7t^7r°, ri' — )► p°y, r\' — )• r|7t+7r~ and w — t- n^n^n^ final states. Pairs of oppositely 
charged particles identified as muons, each having > 550 MeV/c and originating from 
a common vertex, are combined to form J/i[) — t- |J.^|J.' candidates. Well identified muons 
are selected by requiring that the difference in logarithms of the global likelihood of the 
muon hypothesis, Aln£^h; provided by the particle identification detectors 14 , with 



respect to the hadron hypothesis is greater than zero. The fit of the common two- prong 
vertex is required to satisfy < 20. The vertex is deemed to be well separated from 
the reconstructed primary vertex of the pp interaction by requiring the decay length 
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significance to be greater than 3. Finally, the invariant mass of the dimuon combination is 
required to be within ±40 MeV/c^ of the nominal J/i|» mass |15j . 

To identify charged pions the difference between the logarithmic likelihoods of the 
pion and kaon hypotheses provided by RICH detectors, Aln^^K, should be greater than 
zero. In the reconstruction of the B^^^ — )■ J/vJj tt+tt" decay this requirement is tightened 
to be A In C^k > 2 so as to suppress the contamination from B°g^ — )■ J/ij; 7tK decays with 
misidentified kaons. In addition, the pion tracks are required to have px > 250 MeV/ c. A 
minimal value of Axjp, defined as the difference between the of the primary vertex, 
reconstructed with and without the considered track, is required to be larger than four. 

Photons are selected from neutral clusters in the electromagnetic calorimeter with 
minimal transverse energy in excess of 300 MeV. To suppress the large combinatorial 
background from 7r° — )• yy decays, photons that can form part of a 7r° — )• yy candidate with 
invariant mass within ±25 MeV/c^ of the nominal 7t° mass are not used for reconstruction 
of Tj — 7- yy and r|' — )• p°y candidates. 

The ri — t- yy (7t° — ?■ yy) candidates are reconstructed as diphoton combinations with 
invariant mass within ±70 (25) MeV/c^ around the nominal r\ (7t°) mass. To suppress the 
combinatorial background to the r\ — )■ yy decay, the cosine of the decay angle 6'*, between 
the photon momentum in the r\ rest frame and the direction of the Lorentz boost from 
the laboratory frame to the r\ rest frame, is required to have |cos0*| < 0.8. 

The ri' candidates are reconstructed as r\TC^TC and p°y combinations with invariant 
mass within ±60 MeV/c^ from the nominal Tj' mass. For the ri' — )■ p'^y case, the invariant 
mass of the tt+tt" combination is required to be within ±150 MeV/c^ of the p° mass. 
For ri — )■ W^n'n^ (cu — )■ 7i^n~n^) candidates the invariant mass is required to be within 
±50 MeV/c^ of the nominal ri (cu) mass. 

The B°g-| candidates are formed from J/i|;X° pairs with px > 3GeV/c for the X''. To 
improve the invariant mass resolution a kinematic fit [16] is applied. In this fit, constraints 



are applied on the known masses 15 of intermediate resonances, except the wide p" and 
cu states, and it is also required that the candidate's momentum vector points to the 
associated primary vertex. The per degree of freedom for this fit is required to be less 
than five. Finally, the decay time (ct) of the B^^^ candidates is required to be in excess 
of 150 [xm. 

4 Evidence for the J/i() cu decay 

The invariant mass distribution of the selected i/^\> cu candidates is shown in Fig. [2| where 
a B° signal is visible. To determine the signal yield, an unbinned maximum likelihood fit is 
performed to this distribution. The signal is modelled by a Gaussian distribution and the 
background by an exponential function. The peak position is found to be 5284 ±5 MeV/c^, 



which is consistent with the nominal B*^ mass 15 and the resolution is in good agreement 
with the prediction from simulation. The event yield is determined to be y-^o = 72 ± 15. 

The statistical significance for the observed signal is determined as 
S = X \n{CB/ jCs-^-b), where Cs+b and Cb denote the likelihood of the sig- 
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5.1 5.2 5.3 5.4 5.5 5.6 

nij/^^ [GeV/c^] 

Figure 2: Invariant mass distribution for selected — )• J/ij) cu candidates. The black dots 
correspond to the data distribution, the thick solid blue line is the total fit function, the 
blue dashed line shows the background contribution and the orange thin line is the signal 
component of the fit function. 



nal plus background hypothesis and the background hypothesis, respectively. The 
statistical significance of the signal is found to be 5.0 standard deviations. Taking into 
account the systematic uncertainty related to the fit function, which is discussed in detail 
in Section 7.1, the significance is 4.6ct; this also takes into account the freedom in the 
peak position and width in the nominal fit. 

To demonstrate that the signal originates from B° — J/i|;tu decays, the sPlot tech- 
nique [rrj has been applied. Using the J/i[) Ti+Tt'yy invariant mass as the discriminating 
variable, the distributions for the invariant masses of the intermediate resonances 7T° — )■ yy 
and cu — t- 7r+7t^7r° have been obtained. The invariant mass window for each corresponding 
resonance is released and the mass constraint is removed. 

The invariant mass distributions for yy and 7t"'"7t~7r° from B*^ — )■ J/i[) cu candidates 
are shown in Fig. [s} Clear signals are seen for both the cu — )■ 7t+7t^7r° and 7r° — )■ yy 
decays. The yy distribution is described by a sum of a Gaussian function and a constant. 



4 



> 



w 

o 



60 
50 
40 
30 
20 
10 

-10 



_ 1 1 1 1 1 1 , 1 1 1 1 

- 7r° — 7- yy 

/ N 1 / \ 


1 1 1 1 1 _ 

LHCb J 


- , , 1 , , , 1 , , , 1 


1 , , , 1 - 



100 



120 140 



160 180 

[MeV/c^] 




850 900 

[MeV/c2] 



Figure 3: Background-subtracted (a) yy and (b) n^n^T^ invariant mass distributions for 
B° — )■ J/4)7r'*'7r~77 decays. In both distributions the hne is the result of the fit described 
in the text. 



The cu —7- n^n^TfP signal is modelled by a convolution of a Gaussian and a Breit-Wigner 
function with a constant background. The peak positions are in good agreement with the 
nominal tt" and cu masses and the yields determined from the fits are compatible with the 
B° — !■ J/vJj uj yield. The nonresonant contribution in each case is found to be consistent 
with zero. 



5 Decays into J/i()r|^'^ final states 

The invariant mass spectra for Bg — t- J/ij^ri'^'^ candidates are shown in Fig. |4| where 
signals are visible. To determine the signal yields, unbinned maximum likelihood fits are 
performed. For all modes apart from JAprj' (ri' — t- p^y), the B^ signal is modelled by a 
single Gaussian function. In all cases there is a possible corresponding B° signal, which 
is included in the fit model as an additional Gaussian component. The difference of the 
means of the two Gaussians is fixed to the known difference between the Bg and the B° 
masses 18 . Simulation studies for the J/ipTl' ("H' P°y) niode indicate that in this case 
a double Gaussian resolution model is more appropriate. The mean values of the two 
Gaussian functions are required to be the same, and the ratio of their resolutions and the 
fraction of the event yield carried by each of the Gaussian functions are fixed at the values 
obtained from simulation. 

The combinatorial background is modelled by an exponential function. In addition, a 
component is added to describe the contribution from partially reconstructed B decays. It 
is described with the phase space function for two particles in a three body decay under 
the hypothesis of B — )■ J/i|)ri^'''X decay, where X can be either a kaon or a pion, which 
escapes detection. The phase space function is convolved with a resolution factor, which is 
fixed at the value of the signal resolution. 
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Figure 4: Invariant mass distributions for selected Bg 



J/vl^ri^') candidates: 



(a) J/x|;Ti(Ti^yy), (b) B" -> J/tI;ti (t] ^ 71+71-71°), (c) B^ -> J/x^ti' (ti' ^ pV) 

and (d) Bg — )■ J/i[)ri'(ri' — )■ 7t+7t-ri). In all distributions the black dots show the data. 
The thin solid orange lines show the signal B^ contributions and the orange dot-dashed 
lines correspond to the B° contributions. The blue dashed lines show the combinatorial 
background contributions and the dotted blue lines show the partially reconstructed 
background components. The total fit functions are drawn as solid blue lines. The results 
of the fit are described in the text. 



The fit results are summarized in Table [TJ In all cases the position of the signal peak is 
consistent with the nominal B^ mass [Ts] and the resolutions agree with the expectations 
from simulation. The statistical significances of all the B^ decays exceed 7ct. 

To test the resonance structure of the B^ — )■ J/ij^ri'^'^ decays, the sPlot technique is used. 
For the 7t°, r| and r\' candidates the background-subtracted invariant mass distributions are 
studied. The restrictions on the invariant mass for the corresponding resonance are released 
and the mass constraints (if any) removed. The background-subtracted distributions are 
then fitted with the sum of a Gaussian function and a constant component for the resonant 
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Table 1: Signal yields, y^o, the fitted mass, mgo and mass resolutions, a^o for the 
B° ^ J/ipriW decays. 



Mode 3^B0 



niBO ctbo 
[MeV/c2] [MeV/c21 



RO ^ J/i|;ri(ri ^ yy) 810 ± 65 5367.2 ± 3.5 40.1 ±3.6 

B° ^ J/i|;ri(ri ^ 71+71-71°) 94 ± 11 5368.4 ± 2.6 20.3 ± 2.3 

B° ^ J/i|;ri'(ri' ^ p°y) 336 ± 30 5367.0 ± 1.1 8.0 ± 1.1 

RO ^ J/i|;Ti'(ri' ^ 7r+7r--n) 79 ± 10 5369.0 ± 2.8 20.7 ± 2.3 



and nonresonant components respectively. In the fit of the dipion invariant mass for the 
r\' — > Tt+TT^y decay a modified relativistic Breit-Wigner function is used as the signal 



component |19 20 



Background-subtracted invariant mass distributions of the intermediate resonance 
states from the B^ — J/\|;X° decays, are shown in Fig. [5} Clear signals are seen. In all 
cases the signal yields determined from the fits are in agreement with the event yield in 
the Bg signal within one standard deviation (Table [T|. The signal positions are consistent 
with the nominal masses of the r\^'^ mesons and the nonresonant contribution appears to 
be negligible. In each case the invariant mass resolution agrees with the expectation from 
simulation studies. 



6 The J/^\) n+n- decay 



The B° — )■ J/ij'P^ (p° Tt+Tr^) decay is used as a normalization channel 21 . Since it 
contains a J/il^ meson and two pions in the final state, the systematic uncertainty is reduced 
in the ratio of the branching fractions, as the corresponding reconstruction and particle 
identification uncertainties are expected to cancel. 

The invariant mass spectrum for B^^-j — )■ J/i|j7r+7r- candidates is presented in 
Fig. [6| where three clear signals are visible. Two narrow signals correspond to the 
B'^ —7- J/il)7r+7t- and B^ — t- J/i|)7t+7t- decays. The latter decay has been studied in 



detail in Refs. [22, 23 . The peak at lower mass corresponds to contamination from 
B° J/'\\)K*^ (K*° — )■ K~^7i~) decays with a kaon being misreconstructed as a pion. A 
contribution from B^ J/i]; K*° decay is considered to be negligible. 

The invariant mass distribution is fitted with a sum of three Gaussian functions to 
describe the three signals, and an exponential function to represent the background. The 
fit gives a yield of 1143 ± 39 for B° -> J/ip 71+71". 

Previous studies at BaBar 21 show that the B° — )■ J/il;7r+7r~ final state has contribu- 



tions from decays of p° and Kg mesons, as well as a broad S-wave component. A further 
component from the f2(1270) resonance is also hinted at in the BaBar study. To study the 
dipion mass distribution the sPlot technique is used. With the J/4)7t+7t~ invariant mass 
as the discriminating variable, the 7r+7t~ invariant mass spectrum from B° — > J/i|)7r+7t~ 



7 




650 



1000 



650 



Figure 5: Background-subtracted invariant mass distributions for (a) yy from 
Bg — > J/il;ri(ri — )■ yy); (b) Ti+TfTt^ from B^ — > J/i|;ri(r| — )■ 7t+7t~7r°); (c) and (d) 7t+7r~y 
and 7t+7t^ from B^ — )• J/v|)ri'(r|' — )■ p°y, p — )■ tt+tt^); (e) and (f) rjTr+Tr" and yy from 
Bg — 7- J/vJ;ri'(ri' — )• ri7t+7t^). The purple line is the result of the fit described in the text. 



decays is obtained (see Fig. [?]). A dominant p° signal is observed together with a narrow 
peak around 498 MeV/c^ due to Kg decays. There is also a wide enhancement at a mass 
close to 1260 MeV/c^. The position and width of this structure are consistent with the 
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Figure 6: Invariant mass distribution for selected B^^^ — )■ J/i|)7t"'"7t candidates. The black 
dots show the data. The dot-dashed thin orange line shows the signal B° contribution 
and the orange solid line shows the signal contribution, a reflection from misidentified 
B° — )• J/i[) (K* — )■ Ktt) is shown by a blue dotted line. The blue dashed line shows the 
background contribution. The total fit function is shown as a solid blue line. 



interpretation as a contribution from the f2(1270) state. This will be the subject of a 
future publication. 

The distribution is fitted with the sum of several components. A P-wave modified 
relativistic Breit-Wigner function 19 20 multiplied by a phase space factor describes the p° 
signal. A D-wave relativistic Breit-Wigner function is added to describe the enhancement 
at 1260 MeV/c^. The parameters (width and mean value) of this function are fixed to 
the known f2(1270) mass and decay width [T^. The S-wave contribution expected from 
the fo(500) resonance is modelled by a Zou-Bugg 24 ,25 function with parameters from 
Ref. [26]. The p° parameters (mass and width) are fixed at their nominal values and the 
region around the Kg peak is excluded from the fit. The excluded region is ±40 MeV/c^ 
which is four times the mass resolution. A small systematic uncertainty is induced by 
neglecting the p° — cu interference. The value of the uncertainty is estimated to be 0.5% 
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niTt+TT- [GeV/c^] 

Figure 7: Background-subtracted tt+tt^ invariant mass distribution from 
B° —7- J/ipTt+TT^ decays. The black dots show the data. A violet solid line de- 
notes the total fit function, the solid orange line shows the p° signal contribution and the 
blue dashed line shows the f2(1270) contribution. The blue dot-dashed line shows the 
contribution from the fo(500). The region ±40 MeV/c^ around the Kg mass is excluded 
from the fit. 

relative to the p° event yield. 

The reconstruction and selection efficiency for the dipion system has some dependence 
on the dipion invariant mass. A study using simulated data has shown that with the 
increase of the Tf^Tt" invariant mass in the range 300 - 1500 MeV/c^ the efficiency decreases 
by approximately 16%. As the p° meson has a significant width, this dependence needs 
to be accounted for in the determination of the p" signal yield. For this, the efficiency 
dependence on tt+tt" invariant mass extracted from the simulation is described with a 
linear function. Then each entry in the invariant mass distribution is given a weight 
proportional to the inverse value of the efficiency function and the efficiency-corrected 
invariant mass distribution is refitted with the same sum of functions to extract the 
efficiency-corrected event yield for B'^ — > J/ipp*^. The resulting fit parameters both for the 
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Table 2: Fitted yields of the p° resonance, the relative yields of the f2(1270) and fo(500) 
components and probabilities, P, of the fits to the uncorrected and efficiency-corrected 
7r+7t~ invariant mass distributions. 

Uncorrected fit Efficiency-corrected fit 
event yield 811 ±38 (27.6 ± 1.3) x 10^ 

fo (500) fraction 0.20 ± 0.04 0.19 ± 0.04 

£2 (1270) fraction 0.14 ± 0.03 0.16 ± 0.04 

r \%] 40 46 



uncorrected and efficiency-corrected distributions are listed in Table [2j 

7 Measurements of ratios of branching fractions 

Ratios of branching fractions are measured using the formula 

B xo _ g (B ^ J/il^X°) ^ y (B ^ J/r|;X°) ^ ^bXj/^yq 
B-Y» - B{B^ J/rP YO) 3^ (B ^ J/rj; ¥«) ^ B^o e%%^^^, ' 

where y are the measured event yields, £:*°* are the total efficiencies, excluding the branching 
fractions of light mesons and Bxo{Syo) is the relevant branching ratio of the light meson 



X°(Y°) to the final state under consideration 15 . In cases where decays of different types 



of B mesons are compared, the ratio of the branching fractions is multiplied by the ratio 
of the corresponding b-quark hadronization fractions /d//s [27] . 

The total efficiencies consist of three components: the geometrical acceptance of the 
detector, the reconstruction and selection efficiency and the trigger efficiency. For the 
B'^ —7- J/i|) p° decay, the event yield y implies the value weighted by the selection and 
reconstruction efficiency from Table |2] Only the acceptance and trigger efficiencies are 
included in ^BO%j/\|)pO- efficiency components have been determined using the simulation 
and the values are listed in Table |3l 

For channels with photons and neutral pions in the final states, the reconstruction and 
selection efficiencies are corrected for the difference in the photon reconstruction between 
the data and simulation. This correction factor has been determined by comparing the 
relative yields of the reconstructed B+ -> J/v|)K*+(K*+ K+7t°) and B+ J/il;K+ 
decays. The results of these studies are convolved with the background subtracted photon 
momentum spectra to give the correction factor for each channel. The values of the 
correction factors {r]'^°'^^) are also listed in Table |3] 

7.1 Systematic uncertainties 

Most systematic uncertainties cancel in the branching fraction ratios, in particular, those 
related to the muon and J/4> reconstruction and identification. For the final states with 
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Table 3: Branching fractions of the intermediate resonances, total efficiencies (excluding the 
branching fractions of the intermediate resonances), 5*°*, and the photon and tjP efficiency 
correction factors 7]'^°^'^ for various channels. For the B" — )■ J/4) p'^ decay the total efficiency 
includes only the detector acceptance and trigger efficiencies, as the reconstruction and 
selection efficiency for this channel has been discussed in Sect. |6| 



Mode 13 [%] e'°' [%] 7]"°" 



70 



BO -> J/i|;ri (ri ^ yy) 39.31 ± 0.20 0.236 ± 0.006 98.0 ± 7.5 

BO ^ J/i|;ti (ri ^ Tt+Tr-Tr") 22.74 ±0.28 0.059 ± 0.002 94.1 ± 7.5 

BO ^ J/i|;ri' (ri' ^ pOy) 29.3 ±0.6 0.142 ±0.004 98.0 ± 3.7 

BO ^ J/i|;ri' (ri' ^ 7t+7t-ri) 18.6 ±0.3 0.068 ± 0.003 96.0 ± 7.5 

BO ^ J/i|;cu (a) ^ 71+71-71°) 89.2 ±0.7 0.043 ± 0.002 94.1 ± 7.5 

BO ^ J/i|;pO(pO ^ 7t+7r-) 98.90 ±0.16 12.6 ± 0.5 



photons the largest systematic uncertainty is related to the efficiency of 7r0/y reconstruction 
and identification, as described above. The uncertainties of the applied corrections reflect 
simulation statistics, and are taken as systematic uncertainties on the branching fractions 
ratios. 

Another systematic uncertainty is due to the charged particle reconstruction efficiency 
which has been studied through a comparison between data and simulation. For the ratios 
where this does not cancel exactly, the corresponding systematic uncertainty is taken to 



be 1.8% per pion 28 



The systematic uncertainty related to the trigger efficiency has been obtained by 
comparison of the trigger efficiency ratios in data and simulation for the high yield decay 



mode B^ — )• J/ij) with similar kinematics and the same trigger requirements 29 . This 
uncertainty is taken to be 1.1%. 

In the ratios where decays of B mesons of different types are compared (BO or B^), 
knowledge of the hadronization fraction ratio /d//s is required. The measured value of 
has an asymmetric uncertainty of I75 / 



this ratio 
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/o. 

Systematic uncertainties related to the fit model are estimated using a number of alter- 
native models for the description of the invariant mass distributions. For the Bg — )■ J/\[)ri'-''' 
decays the tested alternatives include a fit without the BO component, a fit with the 
means of the Gaussians fixed to the nominal B meson masses, a fit with the width of 
the Gaussians fixed to the expected mass resolutions from simulation and substitution of 
the exponential background hypothesis with first- and second-order polynomials. This 
uncertainty is calculated for the ratios of the event yields. For each alternative fit model 
the ratio of the event yields is calculated and the systematic uncertainty is then determined 
as the maximum deviation of this ratio from the ratio obtained with the baseline model. 

A similar study is performed for the B^ — J/i|; o) channel. As the fit with one Gaussian 
function is the baseline model in this case, here the alternative model is a fit with two 
Gaussian functions (allowing a possible B^ signal). 

In the BO — > J/ip pO case, an alternative model replaces the Zou-Bugg fo(500) term with 
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Table 4: Relative systematic uncertainties for ratios of the branching fractions (7^) for the 
^ J/xpTiW channels [%]. 



Parameter 




-T-)T|'— !>ri7t+7r 


-pT]'-!>T]7t+7T 






Ti'— !>ri7r+7t 


^Icorr 








3.8 


3.9 


3.9 


7t^ reco 


2 X 1.8 


2 X 1.8 




2 X 1.8 






Trigger 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Fit function 


+3.7 
-0.0 


+9.9 
-0.0 


+1.3 
-5.6 


+3.4 
-0.0 


< 0.1 


+0.0 
-2.8 


B{r\,r\',w) 


1.3 


1.7 


2.0 


2.1 


1.8 


2.6 



Table 5: Systematic uncertainties for ratios of the branching fractions (TZ) relative to 
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Parameter TZ 



BO,ti^YY 

B0,pO^7T+7t- 



^^B0,p0^7T+7r- 



^Bg ,ri'— >-Ti7t+7r 
^^B",p0^7T+7r- 



V 



B^'.CU— >-7t+7t 7t' 







BV 



^Icorr 

7t^ reco 
Trigger 
Fit function 



7.6 
2 X 1.^ 
1.1 

+5.1 
-3.7 

0.5 



8.0 



1.1 

+5.0 
-4.3 

1.2 



3.8 



1.1 

fS.O 
-5.7 

2.1 



7.8 



1.1 

+5.0 
-8.7 

1.6 



8.0 



1.1 

+6.4 
-8.8 

0.8 



a Breit-Wigner shape. The mass and width of the broad fo(500) state are not well known. 
The mass measured by various experiments varies in a range between 400 and 1200 MeV/c^ 
and the measured width ranges between 600 and 1000 MeV/c^ 15 . Therefore, the fo(500) 



parameters are varied in this range and the p° yield is determined. Again, the maximum 
deviation from the baseline model is treated as the systematic uncertainty of the fit. 

The uncertainties related to the knowledge of the branching fractions oi r\, r\' , ttP and 
tu decays are taken from Ref. 15 . Other systematic uncertainties, such as those related 



to the selection criteria are negligible. The systematic uncertainties are summarized in 
Tables |4] and [5} The total systematic uncertainties are estimated using a simulation 
technique (see Sect. 7.2). 



7.2 Results 

gO / gO _(/) gO 

The final ratios TZ-QO -q ? ^b" po ^'^'^ ^Bo'po determined using a procedure that combines 
X^-minimization with constraints and simplified simulation. First, the is minimized 

2 _ \ ^ 2 

X / J Xi 5 

i 

where the sum is performed over the six measured event yields for the six different 
modes: B° -> JAl^ri (r] yy), B° J/ipr] (ri 7t+7r"7r°), B° J/i|;ri' (ri' p°y), 

BO ^ {■<]' ^rin+n-), B" ^ J/iptu and B^ ^ J/i|;pO, and xf = In this 
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procedure the foUowing constraints are imposed 

yBO^J/^|)n(Tl^-n') _ 3^BO^J/T|)Ti(ri-5>7t+7r-7tO) 

^BO^J/^|.rl(rl^YY) X S (rj -> yy) £BO^.y^\,^(^^n+^x-n») X i5 (tI ^ 7t+7t-7t0) 

3^BO-s.J/i|jTi'(ri'->pOy) _ 3^BO^J/4)Ti'(ri'-5>Ti7t+7t- ) 



£BO^J/i|)ri'(ri'^pOY) X i3 (ti' ^ pV) £B0^J/i|)ri'(Ti'^ri7t+7t-) X i3 (t]' ^ TlTt+Tt ) 

The ratios T^ro'^ , ^Ro'no' and T^ro'^o are determined using the event yields obtained 
from the minimization procedure. For this determination the efficiencies have been 
varied using a simplified simulation taking into account correlations between the various 
components where appropriate. As both the and the ratios TZ depend only on the 
ratios of efficiencies, systematic uncertainties are minimized. The remaining systematic 
uncertainties have been taken into account as uncertainties in the efficiency ratios. In 
total, 10^ simulated experiments with different settings of Ei have been performed. The 
symmetric 68% intervals have been assigned as the systematic uncertainty. 
The obtained ratios TZ are 

nS,'"^ = 0.90 ±0.09 

B^,T1 

nfo'2o = 3.75 ±0.31 

7^B=o;pO = 3.38 ±0.30 
= 0.89 ±0.19 

where the first uncertainty is statistical and the second is systematic. 
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8 Summary 

With 1.0 fb^^ of data, collected in 2011 with the LHCb detector, the first evidence for the 
B*^ —7- J/i|) cu decay has been found, and its branching fraction, normalized to that of the 
B° —7- J/ip p° decay, is measured to be 

fg;^^^ 0.89 ± 0.19 (=ut)«S(=yBt). 

Multiplying by the known value of i3(B° J/ip p°) = (2.7 ± 0.4) x 10"^ |2l], the absolute 
value of the branching fraction is 



B{B" 3/\\)w) = (2.41 ± 0.52 (stat) Z'^f^ (syst) ± 0.36 (i3B0^j/^,,po)) x 10"^ 

Using the same dataset, the ratio of the branching fractions of Bg — t- JAjTi and 
Bg —7- J/i^ri' decays has been measured. As each of the decays has been reconstructed in 
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two final states, the resulting ratio has been calculated through an averaging procedure to 
be 



0.90 ±0.09 (stat) +|J:[J^ (syst). 



i3(B0 ^ J/i|;ti) 

This result is consistent with the previous Belle measurement of 'R-qV^ = 0.73 ± 0.14 3 , 
but is more precise. Assuming that the contribution from the purely gluonic component is 
negligible, this ratio corresponds to a value of the t] — t]' mixing phase of 4)p = (45.5 ■ 
The branching fractions of the Bg — J/ipr] and B^ — )■ J/i^ri' decays have been determined 
by normalization to the B" — J/ij) p° decay branching fraction, and using the known value 

of /s//d ~ 



0.267 



W.021 
-0.020 



27 



their ratios are 



g(B° ^ J/i1;ti) 

g(B° ^ 

i3(B0 ^ J/i|; pO) 



14.0 ± 1.2 (stat) til (syst) 



1.0 



1.0 



12.7 ± 1.1 (stat) 1° J (syst)_0 9 



When multiplying by the known value of B{B^ — ?■ J/ipp''), the branching fractions are 
measured as 



i3(B° ^ 
i3(B° ^ J/iPV) 



3.79 ± 0.31 (stat) Igf (syst) t'oM (j) ± 0.56 (i^Bo^j/^po; 
3.42 ± 0.30 (stat) +0;^^ (syst) t'olt (f) ± 0-51 [Bbo^j/^.o] 



X 10-^ 
X 10-^ 



The branching fractions measured here correspond to the time integrated quantities, 
while theory predictions usually refer to the branching fractions at t = 0. Special care needs 
to be taken when the Bg and B° decays are compared at the amplitude level, corresponding 
to the branching ratio at t = 30 . Since the J/ipri*^') final states are CP-eigenstates, 



the size of this effect can be as large as 10%, and can be corrected for using input from 
theory or determined from effective lifetime measurements 30 . With a larger dataset such 
measurements, as well as studies of r| — rj' mixing and measurements of CP asymmetries 
in the Bg — ?■ J/ipri*^'^ modes will be possible. 
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